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£Nj ' Abstract 

An influence of nonlinear interference processes at quantum transitions under strong resonance electro- 
magnetic fields on absorption, amplification and refractive indices as well as on four-wave mixing processes 
is investigated. Doppler broadening of the coupled transitions, incoherent excitation, relaxation processes, 
as well as power saturation processes associated with the coupled levels are taken into account. Both closed 
(ground state is involved) and open (only excited states are involved) energy level configurations are con- 
sidered. Common expressions are obtained which allow one to analyze the optical characteristics (including 
gain without inversion and enhanced refractive index at vanishing absorption) for various V, A and H con- 
figurations of interfering transitions by a simple substitution of parameters. Similar expressions for resonant 
four-wave mixing in Raman configurations are derived too. Crucial role of Doppler broadening is shown. 
The theory is applied to numerical analysis of some recent and potential experiments. 

Keywords: quantum interference, lasing without inversion, resonant four-wave mixing, Doppler and 
V/-) ' strong field effects. PACS: 42.50.Gy, 42.55.-f, 42.65 
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,£h ' Many concepts of quantum optics were originated proceeding from the assumed equality of probabilities of in- 

duced transitions accompanied by an absorption and emission of photons predicted by A. Einstein. Requirement 
of population inversion for lasing is direct consequence from this equality. 

At the presence of several resonant electromagnetic fields, probability amplitudes of a coupled quantum 
states contain several oscillating components at close frequencies. Therefore alongside with squared modules of 
appropriate components cross terms indicating an interference of quantum transitions appear while calculating 
transitions probabilities. The coherent nonlinear optical phenomena stipulated by the indicated evolution of 
quantum states, driven by several fields, were called as nonlinear interference effects (NIE) JlL E|. In quantum 
optics NIE may result in different coupling of a radiation with atoms in absorbing and emitting states controlled 
by the auxiliary fields (3,0]. Various appearances of these effects are feasible. Soon after discovery of lasers 
Rautian and Sobelman [ 5| showed feasibility of amplification without inversion (AW I) in two-level systems. 
The features of AW I in optical three- level systems were explored in H ^|. Studies of NIE in absorption/gain 
spectra including experiments on generation of an optical radiation in three-level systems, so that generation 
was possible only at the expense of nonlinear interference effects, drew much attention in 60th and 70th [Q. 
(Review of relevant optical experiments and of earlier papers on AW I in microwave range see in [^|, || .) Coherent 
population trapping (CPT) is one of the appearances of NIE for the states with negligible relaxation rates 
and Doppler effects. In 80th - 90th studies of coherent interference processes at quantum transitions have been 
attracting much interest again in the context of AW I, electromagnetically - induced transparency (EIT), CPT, 
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enhanced nonlinear optical frequency conversion and other manifestations of these effects |10 . 

In classical terms an emission and absorption of a radiation are stipulated by forced oscillations of bound 
charges and depend on phase difference between radiation and induced oscillations. However, a radiation may 
simultaneously drive several coherent interfering oscillations of a various origin. Depending on a relation of 
their phases and amplitudes the interference can be either constructive or destructive, full or partial. Thus the 
matching components of an optical response can either amplify or suppress each other. On the other hand, 
the macroscopic response of a medium can be thought as result of quantum transitions, at which the photons can 

llth International Vavilov Conference on Nonlinear Optics, 24-28 June 1997, Novosibirsk, Russia SPIE Vol. 3485 



252 



simultaneously contribute in several quantum pathways. By applying semi-classical approach a deep analogy 
with many well known effects of classical physics can be used to interpret and foresee the relevant quantum 
optics effects. Thus, leaving aside classifications of involved elementary quantum processes (introduced and 
valid for weak fields in the limits of perturbation theory) , it is possible to predict and to explain wide range of 
optical processes, stipulated by quantum interference, some of them are quite unusual. 

The objective of the paper is to consider various appearance of interference effects in resonance nonlinear - 
optical processes with the aid of the outlined approach in a context of some recent experiments. The amplitude 
and phase relations of interfering intra-atomic oscillations depend on configuration and on relaxation charac- 
teristics of the coupled transitions, on type of nonlinear-optical process, on intensities and frequency detunings 
of the radiations from resonances. Due to the difference in Doppler shifts the contributions to the macroscopic 
polarization of atoms at various velocities in gases may interfere in a different way too. The interference appears 
differently in an absorption, refraction and in different four-wave mixing {FWM) processes. 

As an illustration of NIE the following results will be presented: 

1 . The possibility of an amplification of a radiation without inversion of saturated populations on resonant 
transition is investigated. Influence of the growth of intensity of an amplified radiation on inversionless amplifi- 
cation in various open and closed transition configurations is analyzed. The conditions are formulated and with 
the concrete examples is shown, that by proper change of incoherent excitation rate of levels and of auxiliary 
radiation intensity the index of an inversionless amplification does not decrease with growth of intensity of an 
amplified radiation. The elements of the theory of such lasing without inversion are presented. 

2. It is shown, that due to Maxwell velocity distribution of atoms and corresponding inhomogeneous broad- 
ening of the coupled transitions, incoherent excitation of the intermediate levels may drastically change both 
spectral properties and a magnitude (by orders of magnitude) of nonlinear susceptibilities for resonant FWM 
processes. As the consequence, important power saturation effects appear. These features must be taken into 
account for explanation of the experiments and optimization of frequency-conversion. Resonant FWM coupling 
of two strong and two weak radiations is considered. Formulas for both cases of coupling, one is relevant to 
coherent population trapping, another - when each level is coupled to only one driving field are derived. The 
outcomes are applied to numerical analysis and to discussion of recent experiments |TT) . 

II ABSORPTION AND REFRACTION INDICES FOR A STRONG RADIATION AT THE PRES- 
ENCE OF OTHER STRONG RADIATION, COUPLED TO AN ADJACENT TRANSITION 

j- First, consider interaction of two strong laser fields at the three- 

dj level system. Possible configurations of such systems are shown in 
FIG.|l|: folded V and A, and cascade - H. In further we shall in- 
vestigate spectral features of a gas material for a radiation E4 at 

g frequency W4, tunable in the vicinity of a transition I — m. It's inten- 

sity is not supposed weak. Depending on the configuration of tran- 
1 sitions under consideration one of the auxiliary strong fields E\, E% 

— or E 2 with frequencies u>i, LO3, o-> 2 , resonant to adjacent transitions 

shown in the figure is turned on. All radiations are supposed to 
FIG. 1: Transition configurations. be uniformly polarized CO- or counter-propagating travelling wave: 

Ej(z, t) = Ej exp{— [i(u>jt — kjz)]} + k.c, where kj - can take both positive and negative values, j — 1, 2, 3, 4. 
Incoherent excitation of the levels with Maxwell's velocity distribution, all possible population and coherence 
relaxation channels are accounted for. 

It is necessary to distinguish the open and closed energy-level configurations. In open one (lowest level is 
not ground), the rate of incoherent excitation of the levels by an external source practically does not depend 
on the rate of induced transitions between considered levels. In the closed one (lowest level is ground one), the 
excitation rate for atoms at different levels and velocities depends on the value and velocity distribution of the 
other populations, which are dependent on the intensity of the driving fields. 

A General equations for absorption end refraction indices 

Power dependent susceptibility X4, responsible for absorption and refraction, can be found from the equation: 

P NL (io i )=N Xi E i , (1) 

where polarization P (0J4) is convenient to calculate with aid of density matrix pif. 

P = Np ij d ji + c.c. (2) 



m 

CO, 
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Taking into account above discussed relaxation and incoherent excitation processes, density matrix equations 
for a mixture of pure quantum mechanical ensembles in the interaction representation can be written in general 
form as: 

LnnPnn = <Zn — i[V, p]nn + ImnPmm, ^ImPlra = ^iPi = —i[V, p]lm, (3) 

Lij — d/dt +Tij, Vi m = Gi m ■ exp{i[Q,4t — kz]}, G/ m = — E 4 • di m /2?i, 

where fi 4 = w 4 — u> mn - frequency detuning from resonance; T mn - homogeneous half-widths of transitions, in 
absence of collisions T mn = (T m + T n )/2; T n — ■ 7nj - - inverse lifetimes of levels; 7mn - rate of relaxation from 
the level m to n, q n = Y] j w njTj - rate of incoherent excitation to a state n from underlying levels. For open 
configurations (fo - is mainly determined by the population of the ground state and practically does not depend 
on the driving fields. 

In a steady-state regime a set of density-matrix equation may be reduced to the set of algebraic equations 
|Q. Below we present only results of calculations. Despite of essential distinctions in manifestations of NIE in 
different open and closed configuration, formulas for absorption/gain (a) and resonant part of refractive (Sn) 
indices and also for power dependent populations of the levels can be presented uniformly for all configurations 
shown on FIG.|l|: 



a A /a 0A = Re{xi/x A }, Sn 4 /2Sn Q 4 : = Im{x A /Xi\, a>i/a oi = Re{xi/Xi}, Sn l /25n Qi = Im{x l /x i }, 
XA _ r 4 Ar 4 (l + u 2 ) T Ar t g 2 xj_ _ T l An(l + gp =p Ar 4 tt| 
Xl ~ Pa An 4 (l + 9l + u 2 ) ' X ° ~ Pi An 4 (l + g\ + up ' 



(4) 



Here and further index i specifies transition, resonant to the auxiliary radiation (see FIG .[!]), x i s susceptibility, 
ao,dno,xo are corresponding maximum resonant values at zero field intensities, Pj = Tj + iVLj (for example: 
Pin = Pi = r 4 + iCli, Pim — P m i, Pif — Pi2 = Tif + i(Qi + H2) etc.). If the atom moves with speed v, Doppler 
shift of resonances must be taken into account by substitution flj for = flj — kjv. In further strokes will 
be omitted, but it is supposed, that the Doppler shift in the formulas is taken into account. Ar 4 = r; — r m is 
power dependent population difference; Art 4 = rii — n m , - population of the level in absence of driving fields, 
which is described by the formula: rii = (qi/Ti) + (^fki/^i)(qk/^k)- 

01 = \G t \ 2 /P i P il ,g 2 = \G x \ 2 /P*P 4l , Ul = |G 4 | 2 /P 4 P 4l , U2 = \Gi\ 2 /P*P u . (5) 

Gj are coupling Rabi frequencies. Formulas for populations differences can be presented uniformly too: 

Ar 4 = (An 4 X 2 T An i X 3 )/{X 1 X 2 - X 3 X 4 ), Ar, = {A ni X x T Ati 4 X 4 )/(X 1 X 2 - X 3 X A ). (6) 

r 4 g 2 Ti \ -\- q't 

X 2 = l + i?e{a 24 a5 4 — — ■ + a 2i a? 4 — — — — — -}, 

Pi 1 + 9i + u 2 Pil + gl+u 2 

v u r r4 92 . I + g\ 

X 3 = J Re{a 34 ae 4 — — ■ + a 3i £c t — - — -}, 

Pi 1 + gi + u 2 Pil + g* + u 2 

r 4 1 + u 2 Ti u\ 

Xi = 1 + i?e{ai 4 a? 4 — — ■ + a H ffi 4 — — — — — -}, 

Pi 1 + gi + u 2 Pi 1 + g* + u 2 

r 4 1 + u 2 Ti u\ 

Xi = i?e{a 44 ac 4 — — ■ h ausei — — - — — — -}, 

Pi 1 + g x + u 2 P, 1 + g* + u* 2 

an — — a 34 , an — —a 3i , a 24 = — a 44 , a 2i = —an. (7) 

Sign minus in concerns to folded V (Ei,Ep and A {Ei,E 3 ) schemes, plus - to cascade H {Ei,E 2 ) 

scheme. Beside that in the ladder ii-scheme, Pi must be substituted for P* , and vice a versa: P* for Pj. ae 4 
and c&i - are saturation parameters accordingly for transitions 4 and i. For open configurations 

ae 4 = 2|G 4 | 2 (r ; + v m - 7 4)/(r i r m r 4 ), (8) 

whereas aei and parameters ay depending only on relaxation constants are defined below for each configuration. 

1. V - scheme (fields P 4 , E\; i = 1) 

OPEN CONFIGURATION 

2|Gi| 2 (r g + r ; - 7l ) r g - 7 i r m - 74 

ee l = ee 1 = , a 2i = a 21 = 1, an = 1, a 3i = a 31 = - — — , a 44 = . (9) 

r s r z ri r 3 + r ; - 7i r z + r m - 74 
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CLOSED CONFIGURATION 



ae 4 = 4|G 4 | 2 /r m r 4 ,as l = sei = 4|G 1 | 2 /r g r 1 , 
a 3i = a 3i = O.5A714, 044 = 0.5Ani, a-n = a 2 \ = 0.5[1 + Ani], 044 = 0.5[1 + A714]. (10) 

2. A - scheme (fields E A , E 3 , i = 3) 

OPEN CONFIGURATION 

ac 3 = 2|G 3 | 2 (r m + T„ - 7 3 )/r m r„r 3 , a 2i = a 23 = l, au = 1, 

a-3i = a 33 = r n (r ; - 7 4 )/r/(r m + r„ - 73), a 44 = r ; (r n - 7 3 )/r„(r m + r ; - 74). (11) 

CLOSED CONFIGURATION 

3C4 = a 3l = a 33 = 1 + An 4 - (1 + 2An 4 ) r ^ 73 , a 44 - 0.5[1 - + An 8 (l + ^)], 

a 2l = a 23 = 1 + An 3 (r„ - r m + j 3 )/(T n + T m - 73 ), a u = 0.5[1 + An 4 (l + J 3 /T n )]. (12) 

H-scheme (fields E4,, E 2 ; i = 2) 

OPEN CONFIGURATION 

2|G 2 | 2 (r / +r m - 72 ) , 1 r ; - 7 4 r f - l2 r ; 

862 = ; aw = 1; «2i = a 22 = 1; = 032 = — = ; a44 = . (13) 

1 /I ml 2 1 J 1 TO + 1 / — 72 1 I + 1 m — 74 

CLOSED CONFIGURATION 

ac 4 = 4|G 4 | 2 /r m r 4 , a 3i = a 32 = (1 + 2An 4 )(r/ - 72 )/(r m + T f - 72) - An 4 , a 44 = 0.5(1 - Ana), 

aw = 0.5(1 + An 4 ), a 2l = a 22 = 1 + An 2 (r ro -T f + j 2 )/{T m +T f - 72). (14) 

NIE are associated with coherence at two-photon transitions and disappear at \Pm \ — > 00. At G4 = formulas 
@> © converge into those, similar to discussed and analyzed in || [| |||9). Following Q, range of parameters 
where amplification is not accompanied by the inversion of power-dependent populations can be easily found 
from (Q), (|6|). For the resonant coupling in V and A schemes conditions for AWI at the transitions 4 and i take 
the form, correspondingly: 

Ar 4 /Ar, <g 2 /(l + u 2 ); ArJAn < u x /(l + 91). (15) 

Similar formulas for H schemes show, that on the contrary to the previous configurations, inversion of popula- 
tions on the adjacent transition is required for AWI in center of the resonance, or amplification under certain 
conditions arises in wings of a resonance. (More detail formulas and analysis are given in [Q. 
Threshold and output power of lasing without inversion can be found from the equation: 

«4 = T, (16) 

Where T is loss of a radiation from a laser cavity per one pass, scaled to the unit of length of the amplifying 
medium. Thus, the derived expressions determine conditions and characteristics of inversionless generation too. 



B Numerical analysis 

Below we shall apply the derived expressions for numerical analysis of NIE in open (Ne) and closed (Na) 
configurations of transitions. The same transitions of Ne were considered in Q for illustration of possible AWI 
of weak probe field. The formulas for velocity averaged absorption index were derived and difference in NIE 
for backward and forward waves in inhmogeneously broadened transitions was analyzed in || for the cases 
of weak probe field and coupling Rabi frequency of driving field not exceeding Doppler width of the transition. 
Therefore in further main attention will be given to effects accompanying increase of intensity of an amplifying 
radiation. 

FIG. ||a shows that inhomogeneous broadening does not destroy macroscopic coherence effects and NIE. 
The relative change of absorption index by an auxiliary field appears even larger, than for a homogeneously 
broadened transition. Coherent coupling gives rise to amplification at W4 and establishes populations so that 
there is no inversion neither on one- nor on two-photon transitions FIG. ||b. With the growth of E4 populations 
of to and g levels aim to equal magnitude, which indicates appearance of coherent population trapping. In the 
latter case a small modulation of velocity distribution on the level m appears. Like in ||, analysis shows 
that in order to attain AWI certain ratio between initial (in the absence of the radiations) populations must 
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FIG. 2: Velocity-averaged absorption/gain index and populations vs velocities at inhomogeneously broadened I — m 
transition (I - excited state, open configuration) in the presence of strong field at I — g transition (Ne, 2s2 - 2p4 - 
3s2) (co-propagating waves). Ai = 1150nm, A4 = 630nm (Ne), F m = 3 ■ 10 c _1 , F = 5 • 10 7 e _1 , F g — 10 7 c _1 , 
Imi = 7m! = 0.5 ■ 10 7 c . Ratio of initial populations is: TV; : N g : N m — 100 : 50 : 85. a - absorption index 
(j/4 = Q.4,/Yi m ), b - populations vs velocities (z — v/v, frequency W4 is locked to the absorption minimum). 1 - 
Si = |£i£W2fi| 2 /rir 9m = 5 ; s 4 = \E 4 d ml /2h\ 2 /r 4 r 3m = 0. Or = 350r ig ; 2 - Si = 5, S 4 = 1, Oi = 360r i9 . c - 
velocity-averaged absorption/gain index vs field intensities (fli = Q4 = 0). 



be fulfilled. FIG. |]c shows that absorption (gain) strongly depends on the intensities of both driving and 
probe fields. As it was outlined, the open and closed systems differ both in possible magnitudes of relaxation 
parameters and in dependence of incoherent excitation on an intensity of the driving fields. FIG. || considers 
the case, when 36% of a ground state atoms are initially excited by an incoherent pump to a level m, that 
still may correspond to strong absorption at the transition ml. By that strong driving field E\ may produce 
AWI for co-propagating shorter- wavelength weak radiation at u)^, which makes approximately 50% from an 
initial absorption (FIG. [|a). The amplification happens in absence of saturated population inversion for all 
transitions (FIG. |b). It is essential that a population of a top level m depends on the strength of E\ even at 
zero intensities of a probe radiation £4. 
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-200 2CD y 4 4111 u _o.2 -0.1 0.1 z 0.2 

FIG. 3: Velocity-averaged absorption - gain index (a) 
and populations vs velocities (z = v/v) (b) at inhomo- 
geneously broadened / — m transition in the presence of 
strong field at I — g transition (I - ground state, closed 
configuration, co-propagating waves). A4 = 330nm, 
Ai = 590nm (Na; 3P 1/2 - 3S - 4P 1/2 ). V m = 9 ■ 10 7 c _1 , 
r 9 = 63 • 10 7 c~ 1 , N, : N g : N m = 64 : : 36. a: 1 - 
Si = 10, S 4 = 0, yi = 0; 2 - Si = 10, S 4 = 20, yi = 20. 
b - OJ4 corresponds to: 1 - to the right, 2 - to the left 
point of a zero absorption of the appropriate curve in a . 
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FIG. 4: Velocity-averaged absorption/gain index vs field 
intensities for co- propagating waves (fii = ^4 = 0). a - 
7V ; : N g : N m = 64 : : 36; b - TV, : Ng : N m = 60 : : 
40. Other parameters are like in FIG. H. 



Curve 2 (FIG. ||a) shows, that the AWI strongly depend on intensity of an amplified radiation, that is 
accompanied for the given configuration by noticeable change of populations of levels m and I. (FIG. [|b) 
displays energy and velocity distribution of the atoms corresponding to appearance of transparency at ui±. It is 
interesting to note that the distribution sharply varies with increase of intensity of E\ . 

FIG. ^ shows, that AWI may be maintained in a certain level with growth of intensity of an amplified 
radiation by changing incoherent excitation rate and strength of an auxiliary radiation. 
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Ill INTERFERENCE EFFECTS IN RESONANT FWM AT DOPPLER BROADENED TRANSI- 
TIONS 



The use of resonant FWM in gases for frequency conversion allows one to decrease the required power of fun- 
damental radiations down to the magnitudes characteristic for cw lasers jyj [l3|, [l4| [l5|, |l6). FWM concerns 
to so-called coherent nonlinear - optical processes, depending on phase-mismatch. As it was outline above, at 
resonant coupling, various coherent component stipulated by correlated quantum transitions and giving contri- 
bution to the process of radiation conversion can interfere. Constructive interference gives rise to enhancement 
of appropriate components of nonlinear polarization and destructive on the contrary - to elimination. Studies 
of appearances of quantum interference at FWM continue to attract significant interest [jljj in the context of 
possible use for increase of conversion efficiency. The values of interfering components depend on energy level 
populations, which, in turn, depend on intensity of radiations and on processes of incoherent excitation and 
relaxation. Constructive or destructive character of an interference depends on a relation of phases coherent 
component and, therefore, on detunings from resonance and on type of nonlinear - optical process. In gaseous 
media inhomogeneous Doppler broadening of transitions is characteristic of typical experimental conditions. 
Depending on energy level, value and sign of Doppler shift, contributions of atoms, moving at various speeds, to 
macroscopic nonlinear polarization can both enhance and suppress each other. The above listed effects appear 
in a different way in absorption, refraction and in FWM macroscopic polarization. It turns out that even small 
incoherent excitation of the levels at Doppler broadened transitions may drastically change spectral properties 
and by orders of magnitude value of the nonlinear susceptibility |l8|, g, [l9| . The choice of optimal conditions of 
conversion is essentially determined not only by influence of interference processes on nonlinear susceptibilities, 
but also on indices of an absorption (amplification) and refraction for coupled waves propagating through a 
resonance medium. Velocity selective population transfer and other effects of resonant coupling with strong 
fields give rise to specific power saturation effects in FWM . In pd] experimental features were observed, which 
did not find explanations in framework of before published lowest order perturbation theory. 

This section is devoted to investigation of mutual influence of quantum 
interference, relaxation, incoherent pump of levels, Doppler broadening, ef- 
fects of strong fields on FWM nonlinear polarization, and absorption of 
coupled radiations in the context of experiments Jll]] . We shall consider Ra- 
man like coupling (FIG|) and FWM processes of a type lu<± = io\ — u>2 + <^3 
as well as inverse process 103 = L04. — lo\ -\- ui2- That allows to compare in- 
fluence of the above discussed elementary effects on conversion processes in 
different conditions. The elementary quantum mechanical processes, deter- 
mining FWM, essentially depend on that whether the atom at an energy 
level couples with two strong fields or with one strong and one weak field. 
FIG. 5. Coupled transitions. j n ^ e range of negligibly small change of the strong radiations both due 

to absorption and conversion and assuming exact phase-matching, quantum efficiency of conversion E3 in £4 
7794(^2) = (kz/k^E^z) / Ez\ 2 can be presented in the form: 



CDi 



CD2 



C03 



CO4 



ry g4 = hk 4 \ 2TTNx[ 3) E 1 EZ/(Aa/2) \ 2 exp{-a iZ }(exp{-(Aa/2)z} - l]\ 

,(3) 



(17) 



where N - atomic number density, \a ~ nonlinear susceptibility, Aa = a\ + a 2 + &3 — «4 ~ «3 — 0^4, otj 
— absorption index for corresponding radiation. Quantum efficiency for conversion of £4 in £3 is written in 
symmetrical form by substitution of Xi\ a i f° r X3 > a 3- For small length z « min{(a4 i 3) _1 , (Aa/2) -1 } 
spectral features of conversion is determined only by nonlinear susceptibility and by intensities of the strong 
fields. 



%4, 3 = k 3 k 4 I 2nNx ( ilE 1 E 2 



I 2 z 2 



(18) 



A FWM in two strong and one weak fields at the conditions of maximum coher- 
ence and coherent population trapping 

In this section expressions for nonlinear polarization at frequencies uj^ and uj^ will be presented for cases, when 
the fields E\ and £2 are strong. Their frequencies lj\ and L02 are close accordingly to the transition frequencies 
LUig and Lo ng . Radiations £3 and £4 with frequencies UJ3 and u>i, close to the transition frequencies uj nm and 
LOim - are supposed nonperturbatively weak. NIE in two strong fields may give rise to such population transfer 
between levels I and n that population of the intermediate level g change negligibly. Such behavior is similar 
to CPT (for review see ^0|). In more general sense a term CPT is often applied to the processes whereas 
contribution of NIE in population transfer is crucial. Such type of coupling will be considered in this section. 
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Nonlinear FWM susceptibilities X3 4 are determined by the equations: 

P nl (cj 3 ) = Nxi 3) ElE 2 E 4 , P NL (u A ) = Nx^E.E^E, 



(19) 



Expressions for ^3,4 as well as for absorption (refractive) susceptibilities Xii derived with density matrix in 
similar way as that in Section II, are given by pij| : 



-iK 



X3 
X4 



rf 3 (i + g 2 ) 
-iK 

d 4 (l+9i) 



p* 
Ri 



1 



p 



1.2 



Pi Pi 



12 



Pi 2 

P* 



1 

P32 



R2 


i?4 


P2P12 


*" PaPai 




v R3 


Pul) ' 


P3P32 _ 



(20) 



= TiRi/PiAm, Ri = [(1 + 5 2 )An - u\ Ar 2 ]/(1 + <? 2 * + "2)^2 = [(1 + u 2 )Ar 2 - ,g 3 Ar 1 ]/(l + g 2 + u*), 



9i 



51 



i? 3 = 

IG1I 2 
P 41 P* 

\Gi? 
Pz2(Ia 



Ar 3 (l + gi) - u 3 ii|(l - g 3 ) + qiu 2 Ri 



92 



,92 



|Gx| 2 
Pf 2 P 2 
|G 2 ' 2 



l + 9i 
,.93 = 

",93 



7/4 

ll 



i?4 — 



Ar 4 (l + 92) - .9ii?i(l - 94) + 92.92^2 



|G-' 2 



D* p* 

\Gi? 



,94 



,94 



IGii 2 
P41P4 
2 



', Ul 



, U 2 



|G r ' 2 



,An 



\G2l_ 
AniX 2 



1 

|G 2 | 2 
P12P1 
An 2 A 3 



92+54 



, Ar 2 



|G 2 | 2 
P 32 P 2 * 
An 2 Xi 



■, U4 



|G 2 | 2 ^ 
P32P3 

An 1X4 



P 4 id 3 Pi 2 <*4 PiVa AiX 2 -A 3 A 4 

<i 3 = r 3 + «(n 4 - fii + n 2 ), rf 4 = r 4 + i(n x - n 2 + n 3 ), 

rj = % + Ar 2 (6iaeiP 2 + ^ae 2 F 3 ) - An(&ja3iPi + ^ffi 2 P 4 ), 
Xi = 1 + aiffiiFi - a 2 a; 2 F4, X 2 = 1 + a 3 ae 2 F 3 - a 4 aeiF 2 , 
X 3 = a 2 ae 2 F 3 - ai38iP 2 , X4 = a 4 aeiPi - a 3 83 2 F 4 . 



Xi A" 2 — X3X4 



Pi = Re{ 



9 2 



Pi 1 + g* 2 + u 2 



Here and further Ari 



n 



}, P 2 - Pe{ 



,Ar 2 



hi 



Pi l+5 2 * 
n -r g , Ar 3 



— },P 3 = ^{^ 1+ 1 + " 2 M 1,^ = Re{ 
U 2 r 2 1 + <? 2 + W 2 



53 



P 2 * 1 + <? 2 * + 712 



(21) 

}■ 



r m , Ar4 = n — r m , Ari4 = — n m , etc., 77 and n,- are 



population of levels accordingly power dependent and in absence of the fields; aei, ae 2 - saturation parameters; 

&ki - coefficients, determined by relaxation properties of the transitions, which are different for open and 

closed transition configurations and will be defined below; K is constant. 

OPEN CONFIGURATIONS 

For the open system the parameters in (22)are: 



aei = 

r, 



2|Gi| 2 (r 9 + r i - 7 i) 



«2 



r 9 r ; ri 

- 71 



Ti r„ + r„ 



72 



b{ = - 



6™ = -- 



r ( + r 9 

Tg -72 



7i ' 



04 



14 



33 2 

r, r r; 



2|G 2 | 2 (r g + r n - 72 ) 



r 9 r„r 2 



72 



r„ r ; 4 
r„ 



71 
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a\ = a 3 



J- n 



7i 



r„ - 72 



r* + r 3 - 71 
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CLOSED CONFIGURATION 

The corresponding parameters in (B2) take the values: 



aei = 4|Gi| 2 /r g ri, ai = 0.5[1 + Am(l + j 2 /T n )],a 2 = 1 + Am - (1 + 2Am)(r 9 - j 2 )/(T g +T n - 72 ), 



a 3 = l + An 2 [l-2(r 9 -7 2 )/(r fl +r n - 72 )],a 4 



0.5[1 - ( 72 /r„) + An 2 [l + ( 72 /r„)], 
b[/m = b™/n m = (r„ + 7 2)/2r„, b\j ni = b l 2 / ni - -(r„ - r s + 72 )/(r 9 + r„ - 72 ), 

6? = 0.5n„ - (1 - n„) 72 /2r„, b% = (2n n - l)(T fl - 7 2 )/(r s + T n - 72 ) - n n , 

n„ 



bf = n g (T n + 72 )/2r„ - 0.5, b 9 2 = 1 - n g - (1 - 2n 9 )(r s - j 2 )/{T g + T n - 72 ), 
ni = (1 + w m /T m + w g /r g +w' n /T n y ,n m = w m ni/T m ,n g = w g ni/T g , 
n n = w' n ni/T n , w' n = w n + w 3 7 3 „/r 9 + w m 7„ m /r m 
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B FWM of two strong and two weak radiations under condition of perturbation 
of each resonant energy level only by one strong radiation 

In jllj] the features coming out from increase of intensity of a radiation at frequency to 3 were investigated too. 
Consider a case, where the fields E\ and E 3 are strong, and E 2 and E4 - weak. With the aid of solution of a 
set of equations for off- and diagonal elements of density matrix up to the first order of perturbation theory in 
respect of the weak fields equations for the susceptibilities can be presented as |23] : 



-IK 



X2 = 
X4 
i?2 = 



<fc(l 



v 5 H 

-iK 



3 5 *) 



An Ar 3 



Ar 2 {\ + 37 + v 7 ) 



PiPii 

An 
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Ar 3 



p: \ p, 
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1 
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43 
1 
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,V 5 = 



|G 3 | 2 IGil 2 |G 3 | 2 



P\P\2 P3P32 

V3O- + V7 - gs)Ar 3 - g 3 (l + g 7 - v$)Ari 

(1 + g 2 + v 2 ) + [37 + 32(37 - v$) + v 7 + v 2 (v 7 - g s )] 



,V7 



Pad* 2 

\G 3 \ 2 
P* 2 d\ 



• 95 = 



97 



|G 3 | 2 

P* H* 
r 3 2 u 4 



Ar 4 (l + v 5 + g 5 ) - gi{\ + 35 - v 6 )Ari - n(l + v 5 - g 6 )Ar 3 
Pi = 7— ; s — ; — f ; ; — ; ; n > v 6 



,g& 



(1 + 34 + ^4) + [v5 +V4(v5 - ge) +g5 + 34(35 - we)] ' "° ^43^2' 96 P^id^' 4 P43P4' 
d 2 = Y ng + + fl 3 - fi 4 ), x*/x° = TiAn/PiAm, (i = 1, 3), = T^/pAn^ (i = 2, 4). (23) 

The rest notations are the same as in Subsection III. A. Expressions for the populations are: 
OPEN CONFIGURATION 

An = [(1 + as 3 )Ani + 6ias 3 An 3 ]/[(l + sei)(l + ae 3 ) - aiaei&iae 3 ], 
Ar 3 = [(1 + a3i)An 3 + 0i8eiAm]/[(l + asi)(l + ae 3 ) - aiaei&i8e 3 ], 
Ar 2 — An 2 — 6 2 ae 3 Ar 3 — a 2 33iAn, Ar 4 = An 4 — a 3 ae 1 Ar 1 — 6 3 as 3 Ar 3 , 



P i3 d* 2 
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P* 2 dl 

JAP 
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CLOSED CONFIGURATION 

The populations of levels are described by the equations: 

T m r m = w m r t - 2Re {iG 3 r 3 } , Y g r g = w g n - 2Re {iG*n} , 
r„r„ = w n ri + 2Re {iG 3 r 3 } + 7 s „r g + 7 m „r m , r ; = 1 - r m - r 9 - r„, 

where ri = iG\Ar\/ P\, r 3 = iG 3 Ar 3 / P 3 . The solution is 

n = nj(l + ae 3 )(l + asi)//3, r g = (1 + as 3 )[n;(l + asi) - Ani]/j9, 
r„ = K(l + aes) (1 + asi) + [An 3 (l + asx) + An l72 ae 1 /r„](l + 6a3 3 )} /(3 
r m = {n TO (l + ae 3 )(l + aci) + [An 3 (l + aei) + Arii7 2 a5i/r„]6ae 3 } //3, 
An =n-r g = Am(l + a5 3 )/^, Ar 3 = r„ - r m = [An 3 (l + aei) + An 1 j 2 ee 1 /T n }/ f3, 



ae 3 )[l - An 3 + 2(ni - 

^ m — W>lW rn /T^ r n 1 Tig 



,)aei] + (1 + 26ae 3 )[An 3 (l + 33!) + Ani7 2 a3 1 /r„], Am = n ; 
n/w g /r g , n„ = niWn'/T n , ni = (1 + w m /T m + w g /T g + w n ' ' /T n )~ 



where /3 = (1 
An 3 = n n — rt 
w„' = 

(2|G 3 | 2 (r m + T„ - 7 3 )/r m r n r 3 )(r§/|P 3 | 2 ). The remaining denotations are former. 



jlgn/Tn + W m J mn /T n , b = T n /(T m + T n - 73 ), 35i - (2 1 Gi | 2 /T^g) {Tl/ \P X 1 2 ) , 35 3 = 



C Effect of Doppler broadening on resonant FWM 

Formula for m lowest order of perturbation theory can be derived from (20), (p3) at G 



xf (w 4 = Wi - o; 2 + w 3 ) = 



iii" 



+ i(0 3 -r! 2 ) 



r m / + ^fij - o 2 + n 3 ) 

1 r n g - n n 

~\~ i(fl^ ^ 2 ) ^ n 9 i^ 2 



(24) 
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where f2j = fl, — kj-v, m = JVj • exp{— (v/v) 2 }/t/itv. As the function of v all terms in (p4|), besides those 
proportional to n g — n„, have all poles in one and the same complex half plane. Therefore at I\ << kiV only 
terms, proportional to n g — n n , do not vanish after averaging over Maxwell's velocity distribution. Velocity 
averaged susceptibility is: 

< AS), ) > _ i^ 2 exp{-(n 2 /M 2 }(^ - N n ) 

" feetfi + ^ni-fcina/feJltfa + ^na-fca^/fe)]' 

f i = T n i + (k 1 /k 2 - l)r„ 9 , f 3 = T gm + (fc 3 /fc 2 - l)T ng . 

As it is seen from (^5|), for the process W4 = Wi — lu 2 + ^3 interference of contributions of atoms at different 
velocities to the velocity averaged nonperturbed FWM nonlinear susceptibility < X4 >u leads to the fact, 
that in the lowest order on the small parameter T 2 /k 2 v, it 1S proportional to the velocity integrated difference 
between populations of the excited states N g — N n . In the similar way, with aid of (]2^) one can find, that on the 
contrary, < \2 >v for the process oj 2 = — ^4 +0J3 in the same approximation is determined by the population 
difference on transitions from the lowest level. For the resonant sum frequency FWM 0J4 = uj\ + uj 2 + 0J3 in the 
cascade configuration of levels velocity averaged susceptibility occurs proportional to higher order of the small 
parameter T/kv compared to Raman-type difference- frequency coupling ||. These features demonstrate great 
difference between resonant FWM processes in homogeneously and inhomogeneously broadened transitions. 

In strong electromagnetic fields above mentioned processes are accompanied by the velocity selective pop- 
ulation transfer and by some other intensity dependent effects. In [ pd| experiments on resonant cw FWM at 
Raman-like electronic molecular transitions of Na 2 have been carried out. Frequency tunable radiation at w 2 
was generated at the same transition of Na 2 either in external dimer Raman laser or FWM was performed 
inside the Raman laser cavity. Frequency lo 2 was tuned by tuning uj\. Radiation at L03 was provided from cw 
dye laser. High conversion efficiency have been attained in single frequency nearly power saturation regime. 
Observed FWM frequency tuning characteristics occured in disagreement with the predictions of lowest order 
perturbative theory. From that the authors derived the questions to be answered with the aid of an advanced 
nonperturbative theory. We shall use above presented expressions for the numerical analysis of the models with 
the parameters, close to that in the experiments, in order to explain main observed features. 

The electronic - vibration-rotation transitions between X, A and B electronic levels of the dimer were used 
in the experiments, the lowest electronic level being ground one. Two FWM processes were investigated: when 
frequency u> 3 less than uj 2 and, therefore frequency of a generated radiation W4 was less than w\ (down conversion) 
and opposite upconversion process. As it was discussed above, different appearances of interference processes 
at Doppler broadened transitions can be expected in those cases. Main observed experimental dependencies, 
which did not find explanations, can be summarized as follows. 

According to (|25|), at fii — k\^l 2 jk 2 — 0, in the lowest order of perturbation theory the maximum output of 
FWM at UI4 as a function of UJ3 corresponds to = k 3 ft 2 /k 2 = 0. Lincshape of the resonance is Lorentzian 
with the linewidth of the order of characteristic homogeneous widths of optical transitions. However, in the 
down conversion experiments the wide resonance of the order of Doppler width of transition ml with the center 
being locked at oj 3 « u> mg was observed. It's position practically did not vary at tuning u 2 within Doppler 
resonance of the transition gn (at the expense of tuning of lo\, so that fl 2 — k 2 Vti/ki). In the upconversion 
experiments the resonance was tunable by tuning frequency oj 2 , but with the slope less than dVl^/dQ,^ = k^/k 2 . 
Width of the resonance also was commensurable with the Doppler width of the transition ml. 

For numerical analysis we have used a model with the transitions parameters, close to those from the 
experiment. 

1. Down conversion: A m ; = 598 nm, X g i — 488 nm, A„ m = 655 nm, X gn = 525 nm; kiV = 6.94, k 2 v = 6.45, 
k 3 v = 5.17, k iU = 5.66, (in terms of I0 9 cr 1 ); T m = 200, r„ = 30, T g = 260, j mn = 2, lml = 4, lgn = 20, 
lg i = 10, T ; „ = 40, T nm = 110, Tim = 110, r sm = 130, r ng = 140, r ;3 = 140, (in terms of 10 6 c' 1 ) 
Ni/N n = 30/2. 

2. Up-conversion: A m ; = 473 nm, X g i = 661 nm, X mn = 514 nm, X gn = 746 nm; k\v = 5.12, k 2 v = 4.54, 
k 3 v = 6.59, hv = 7.16, (in terms of 10 9 cr 1 ); T rn = 200, T n = 30, T g = 260, 7mn = 2, 7m/ = 4, lgn = 20, 
lg i = 10, T ln = 40, r„ m = 110, T lm = 110, r sm = 130, T ng = U0, T lg = U0, (in terms of 10 6 cr 1 ), 
Ni/N n = 110/3.2, population of two upper levels being negligibly small. 

First, with an aid of these models we shall illustrate a role of an interference at velocity averaging of nonlinear 
susceptibilities in weak fields. For down conversion in exact one- and multiphoton resonances and homogeneously 
broadened transitions computing gives the ratio of squared modulus of nonlinear susceptibilities |xi /xi | 2 = 
2.5. For averaged values it yields | < ^3 >v / < Xi ^ >v | 2 = 2.31 • 10 2 . If to change the population ratio for 
the inverse magnitude (N n /Ni = 30/2), we obtain: Ixjj^/xi | 2 = 0-4- The difference between averaged values 
sharply decreases. Their ratio in this case yields: | < x 3 3 ' > v / < xi^ >v | 2 = 0.13. 
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For up-conversion similar computations give: \x 3 /X4 \ 2 = 2-7, \ < X3 >v / < X4 >v \ 2 = 1-45 • 10 3 . At 
the inverse population ratio (N n /Ni = 110/3.2) we obtain: \Xa /Xi | 2 = 0.37. The difference between averaged 
values sharply decreases. Their ratio in this case is: | < x^ >v / < X^ >v | 2 = 0.24. 

Thus, effect of inhomogeneous broadening of the resonant transitions on FWM processes may be very 
strongly dependent on a specific process, as well as on distribution of the populations over levels and velocities. 
Therefore one can expect that velocity selective population transfer and other effects of strong fields may change 
conclusions of the lowest order perturbative theory. 

For small number density and medium length FWM conversion efficiency of weak radiation is proportional to 
a product of intensities of the strong radiations and squared modulus of velocity averaged nonlinear susceptibility 
(equation (|l8|)). The later is intensity dependent too. Further we shall numerically analyze effects of the strong 
fields on conversion efficiency with an aid of the expressions of the subsections III. A and III.B. Intensities of 
the radiations will be characterized by the parameters Si = |Gi| 2 /r g ;r n ; and S2 = \G2\ 2 gn Y n i, which are 
chosen in near saturation range like in the experiments. Figures @ show, that for the chosen parameters due 
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FIG. 6: Squared modulus of scaled velocity averaged nonlinear susceptibility | < xi 3 ' >v / < xlo* >« 2 vs scaled 
detuning y 3 = £l'i/T nm at various detunings Hi (il 2 = ta/feiOi. a and b: upconversion, Sj = 150, S2 = 350); a: 1 
- fii = 0, 2 - ill = 10r is , 3 - ill = 20r i9 , 4 - Oi = 301%; b: 1 -fii = 4or i9 , 2 - Hi = 5or i9 , 3 - ill = 60T lg , 4 - 
ill = 70r i9 . c: downconversion; Si = 1000, S 2 = 2000; 1 - ill = 0, 2 - il x = 501%. 
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FIG. 7: Squared modulus of scaled velocity averaged nonlinear susceptibility | < xi 3 ' >v / < xio' >•» | 2 vs scaled 
detuning j/3 = ils/Fnm at various detunings Hi (Q2 = fe/fciHi) (upconversion, each level is coupled to only one strong 
field), a: Oi = 0, Si = 1000; 1 - S 3 = 1000, 2 - S 3 = 1600. b: Si = 2000, S3 = 2100; 1 - yi = 0, 2 - yi = 20, 3 - 
yi = 40, 4 - yi = 60. c: same as in b, but Si = 4000, S3 = 4200. 

to power broadening the resonance is much broader than homogeneous transition width and is commensurable 
with Dopplcr lincwidth, which is of the order 60 in the used scale. For the parameters, corresponding to the 
up-conversion experiments (FIG. ^| a,b) in the range of small detunings Hi (Sl 2 = (FIG. || a) the 

peak of the tuning curve is displaced very insignificantly (and even in the opposite side, depending on the value 
of Hi). At further increase of Hi (FIG. ||b) the maximum shifts with the increase of Hi, so that the slope ^3/^1 
is variable. A maximum of the slope corresponds to the detunings Ox of about a half of Doppler width of the 
transition gl. Thus for the considered intensities the value of the slope reaches « 0.8, that makes w 0. 5(k%/k-i). 

FIG. [| c is computed and drown for the parameters, corresponding to down conversion experiments. For 
the considered intensities the peak occurs locked to the center of transition ml practically in all an interval of 
Hi within the Doppler width of transition gl. 

When the weak field detunings are fixed and driving fields frequencies u>i and LO2 (^2 = (^2/^1)^1) are tuned 
to the maximum, computer analysis of the slope <if2i / dil^ in the range of Si — 65, S2 = 2.33 • Si shows behavior 
similar to that observed in experiments Jll| (dili/dil 3 ps ki/2k 3 ). (Ratio S1/S2 = 2.33 is chosen according to 
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the ratio of prodacts of experimental field intensities and Franck-Kondon factors.) 

In the experiments jll| the features, following increase of intensity of £3, were observed too. In order to 
consider effects of this field and to understand whether CPT play a decisive role in observed dependencies we 
have carried out numerical analysis of the up-conversion model with aid of formulas from the Subsection III.B 
(FIG. §. 

FIG. [7] a shows that at certain ratio of intensities even power narrowing of Doppler broadened FWM 
resonance may happen. FIG. [?] b displays approximately constant slope d£t^/d£L\ (O3 corresponds to the 
maximum output), which is about 0.75, in a quite wide interval of Only for f2i, larger than Doppler width, 
the slope starts decreasing. At larger intensities the slope varies more considerable, when fii is tuned within 
the Doppler line. So for FIG. c the slope makes up 0.1 in the vicinity of y\ = 20; 0.8 - in the range of y\ = 40 
and 0.6 - at y\ = 60. 
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FIG. 8: Velocity-averaged absorption/gain indices scaled to the nonperturbed absorption index at u) — u g i < on > v 
I < a® > v at Si — 150, S2 = 350 (upconversion) . a - absorption of strong radiations vs wi (Q.2 = fe/fcifii); change in 
absorption of converting and generated radiations as tuning u>i (H3 = 0) (b ) and W3 (fii =0) (c ). 

In an absorbing medium spectral properties of absorption indices may 
bring important effect on FWM. FIG. || a-c display corresponding 
lineshapes. FIG. [| is computed with aid of formula ( |l7|) and shows 
that additional broadening of the tuning curve of FWM output may 
arise from the propagation effect in an absorbing medium [Z = aio-z, 
fii = 0, Si = 150, S2 = 350 (upconversion)). 

In conclusion, the theory of nonlinear interference processes at 
Doppler broadened quantum transitions in two strong resonant optical 
fields is developed. The derived formulas allow one account for such 
contributing processes relevant to experiments as various relaxation 
channels, incoherent excitation by an external source, population transfer 
and other coherent and incoherent effects accompanying coupling with 
strong optical fields in various open and closed configuration of tran- 
sitions. Explicit formulae accounting for the effects of the strong fields 
FIG. 9: Effect absorption on FWM. are derived. Such appearance of quantum interference as amplification 
(and lasing) without inversion of power saturated populations and specific effects in resonant four wave mixing 
in gases are analyzed with the aid of numerical models. Parameters of the model are close to those of some 
recently carried out experiments. Crucial effect of Doppler-broadening on the contributions of the populations 
of different levels to four-wave mixing, which determine selection of optimal energy-level configuration and 
conditions for the experiments, is shown. Unexpected dependencies, observed in the experiments are explained. 
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